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^v^j i Abstract. In this paper we consider two a priori very different problems: construction of 

the eigenstates of the spin chains with non parallel boundary magnetic fields and compu- 
tation of the partition function for the trigonometric solid-on-solid (SOS) model with one 
reflecting end and domain wall boundary conditions. We show that these two problems are 
related through a gauge transformation (so-called vertex-face transformation) and can be 
solved using the same dynamical reflection algebras. 
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! 1 Introduction 

O . 

In 1988 E. Sklyanin proposed a way [25] to generalize the algebraic Bethe ansatz [7] to the 
open integrable systems. In particular, this approach permitted to construct the eigenstates 
for the quantum XXZ spin chain with external boundary magnetic fields parallel to the z axis 
(diagonal boundary terms). General algebraic framework of this method (quantum inverse 
scattering method) worked also for the non-parallel boundary magnetic fields, i.e. the double 
^ ■ row monodromy matrix, commuting transfer matrices and the trace identities were obtained in 

the most general case. However, the conservation of the z component of the total spin turned 
out to be essential for the eigenstates construction. 

It is important to mention that the algebraic Bethe ansatz establishes a clear relation between 
the quantum spin chains and two-dimensional models in statistical mechanics. The periodic XXZ 
spin chain was solved using the transfer matrix of the six-vertex model, while the open chain 
with diagonal boundary terms was solved using the transfer matrix of the six-vertex model with 
reflecting ends. 

This equivalence with two-dimensional models in statistical mechanics also turned out to be 
essential for the computation of the correlation functions. It was shown by V. Korepin [20] that 
the partition function of the six- vertex model with domain wall boundary conditions (DWBC) is 
the key element for the study of the correlation functions. The determinant representation for 
this partition function found by A. Izergin [T3] is crucial for the computation of the correlation 
functions starting from the algebraic Bethe ansatz [151 \TE[ [19] . 



*This paper is a contribution to the Proceedings of the International Workshop "Recent Advances in Quantum 
Integrable Systems". The full collection is available at |http://www. emis.de/journals/SIGMA/R AQIS2010. html 
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For the open case the corresponding partition function was computed in [27]. This represen- 
tation was used first to compute the scalar products and norms of the Bethe vectors [25] and 
then to study the correlation functions of the open spin chains with external boundary magnetic 
fields parallel to the z axis [T6l [PTj. 

As we mentioned above, the construction of the eigenstates for the spin chains with arbitrary 
boundary magnetic fields is more complicated. Recent advances in the study of these systems 
permitted to apply the usual Bethe ansatz technique to the non-diagonal case [2TJ O [30] . This 
case has many interesting applications, for example it can be used to study out-of-equilibrium 
systems such as asymmetric simple exclusion processes (ASEP). The Bethe ansatz solution 
permitted to deal with problems which were out of reach for the usual techniques, such as the 
relaxation dynamics [51 [6] . 

The algebraic Bethe ansatz technique developed in [31 [3D] is based on a gauge transformation 
(so-called vertex-face transformation) which diagonalizes the boundary matrices and relates the 
spin chain to a trigonometric solid-on-solid (SOS) model with reflecting ends. Once again, it can 
be demonstrated that these two models can be described using the same algebraic structures. 

The dynamical Yang-Baxter algebra proved to be the main tool to solve the SOS models 
[TOl [9]. In particular, it was used to compute the partition function of the SOS model with 
DWBC [Ml [23]. Unfortunately, the final results obtained in these papers cannot be written as 
a single determinant as in the six vertex case and thus, this result cannot be easily used for the 
computation of the correlation functions. 

In this paper we consider one possible connection between the spin chains with non-diagonal 
boundary terms and trigonometric SOS models in a systematic way. In the first part of our paper 
we review the construction of the eigenstates and we give a more simple picture than [31130] using 
two types of SOS reflection algebras. Finally, we consider the SOS model with one reflecting 
end and we use this algebraic framework to compute the partition functions for different types of 
domain wall boundary conditions generalizing the results of our paper [TT]. The main advantage 
of our result is that all these partition functions can be written as a single determinant. 

The main aim of this paper is to relate the trigonometric SOS model with reflecting end and 
spin chains with non-parallel boundary magnetic fields in a simple and comprehensible way and 
to show the algebraic structures which can be used in the future to compute the correlation 
functions and form factors. 

It is important to note that the Bethe ansatz cannot be applied for arbitrary boundary 
fields. It works only if the parameters satisfy certain boundary constrains. It seems that it is 
absolutely necessary to have at least one relation between the boundary parameters in order to 
construct the Bethe ansatz. Here we consider a special case when two boundary matrices can be 
diagonalized simultaneously by the vertex-face transformation. It requires two constraints [30] 
instead of the only one obtained initially [21] . We would like to underline that this case is more 
complicated than the general one, as in this situation the Bethe ansatz solution is not complete 
(it is a degenerated situation where completeness conjecture [22] does not work). 

The paper is organized as follows. In Section [2] we review the quantum inverse scattering 
method for the spin chains with non-diagonal boundary terms. Following [25] we construct 
the "vertex type" double row monodromy matrix as well as a family of commuting transfer 
matrices. We show how the Hamiltonian can be obtained from these transfer matrices. In Sec- 
tion [3] we introduce the vertex-face transformation which diagonalizes the boundary matrices 
and transform the vertex type monodromy matrices into "SOS type" objects, satisfying dynam- 
ical Yang-Baxter algebra relations. We construct two double row monodromy matrices which 
satisfy the dynamical reflection relations or dual dynamical reflection relations. In Section [5] 
we use these algebras to construct eigenstates of the spin chain with non diagonal boundary 
terms. Finally, in the last section we use these two reflection algebras to compute the partition 
functions of the SOS model with reflecting end and domain wall boundary conditions. Some 
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details on the algebraic Bethe ansatz, Bethe states and relations between the two dynamical 
reflection algebras are given in the Appendices. 

2 XXZ spin chain with non diagonal boundary 

We consider the open XXZ spin chain with the most general non-diagonal boundary terms 

N-l 



H = J2 ( a f a i+i + a i a i+i + c ° sh wf^f+o 



i=i 



H — — T — cosh Q cosh ba\ + sinh raf — i cosh toV\ 

smh sinn o 

smh Tj 

: [ — cosh ( cosh 5a z N + sinh ra% — i cosh Ta^] . (2.1) 



sinh (" sinh 5 



The Hamiltonian acts in the Hilbert space 'K = (8^=1 ^m, ~ C • Two boundary magnetic 
fields are described here by six parameters 5, 5, £, C> t and r. 

This Hamiltonian can be obtained using the Quantum Inverse Scattering Method (QISM) [7] 
or more precisely its boundary version [25] starting from the Yang-Baxter equation and usual 
six-vertex trigonometric i?-matrix. 

The key element of the QISM is the quantum i?-matrix R : C — > End(y <8> V), V ~ C 2 
satisfying the Yang-Baxter equation. 

Rl2{M — M)Rl3{^l — A3)i?23(^2 — A3) = i?23(A2 — Aa)i?i3(Al — X^Rtfi^l ~ X 2 )- 

The solution corresponding to the XXZ spin chain is the six- vertex i?-matrix 



R(X) 



/sinh(A + r?) \ 

sinh A sinh 77 

sinh 77 sinh A 

\ sinh(A + 7/)/ 



This i?-matrix is symmetric, namely VR\ 2 V = R12, where V is the permutation operator on 
V <g> V and satisfies the following important relations: 

• initial condition 

R(0) = sinh?]^, 

• unitarity 

i?i 2 (A)i?2i(-A) = - sinh(A - rj) sinh(A + rj) Id, (2.2) 

• Z2 symmetry 

a\a\R 12 {X)a\al = R 12 (X), 

• crossing symmetry 

- a yR%(-X-r l )a y 1 =R 21 (X). (2.3) 
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Here R 2 \ = V\ 2 R\ 2 V\ 2 , ll denotes the matrix transposition on the first space of the tensor 
product, and <J x ,y,z are the usual Pauli matrices. 

The bulk monodromy matrix To (A) £ End(Vo <8> IK) of the inhomogeneous system is con- 
structed as: 

T (A) = i? i(A - 6) • • • Ron(X - £n)- (2-4) 

In these last expressions, i?om denotes the i?-matrix in End(Vb <S> "Km), and £i,£ 2 , • • • ,€n are 
arbitrary complex parameters (inhomogeneity parameters) attached to the different sites of the 
chain of length N. 

It is easy to show that it satisfies the Yang-Baxter algebra relations: 

fl 12 (Ai - A 2 )T 1 (A 1 )T 2 (A 2 ) = r 2 (A 2 )T 1 (A 1 )i? 12 (A 1 - A 2 ). 

This monodromy matrix can be used to construct the eigenstates of the periodic XXZ 
chain [7]. For the open chain a more complicated construction is needed to reconstruct the 
Hamiltonian and, eventually, to solve the model [25j . 

Taking account boundaries require the boundary K± matrix which satisfy the reflection 
equation [H [25] : 

J R 12 (Ai 2 )iT_(A 1 ) J R 21 (A 12 )^_(A 2 ) = K_(A 2 )iii 2 (A 12 )A'_(A 1 ) J R 21 (A 12 ) (2.5) 
and its dual: 

i?i 2 (A 2 i)^(Ai)i?2i(-Ai 2 - 2 V )K$(\ 2 ) = K%{\ 2 )R 12 {-\ 12 - 2rj)K%(\ 1 )R 21 {\ 21 ), 

here we used the following short notations Ai 2 = Ai — A 2 , Ai 2 = Ai + A 2 . 
We consider the most general solution of these equations [13] : 



( cosh(6 + C)e~ A - cosh(5 - C)e A _ T sinh(2A) \ 

2sinh(5 + A)sinh(A + C) & 2 sinh(5 + A) sinh(A + Q 

r sinh(2A) cosh(<5 + C)e A - cosh(5 - £)e~ A 

V ~ 6 2 sinh(<5 + A) sinh(A + C) 2 sinh(<5 + A) sinh(A + () / 

and K + (X) = K—{— A — 77; 5, C,r). The complete inhomogeneous spin chain algebra is actually 
described by the double row monodromy matrix: 



(C/_)o(A)=r (A)(K„)o(A)r (A) 



or, equivalently, by the dual one: 



A_(A) -B-(A) 
C_(A) D_(A) 



(^)o(A) = T*"(A)(At)o(A)f ^(A) = (j^j 



where, due to the unitarity (12.21) and the crossing symmetry p. 31) : 

f (A) = R m (X + Cjv) • • • Rio(X + £1) = 7(A)^T '°(-A - rj)a y = %\)T~\-\), 
with normalization coefficients: 

N 

7 (A) = (-1)* 7 (A) = 7(A) J] sinh ( A + & " ? ?) sinh ( A + 6 + ^ 

i=l 
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The double row monodromy matrix satisfies the reflection algebra relation 

i?i 2 (Ai 2 )(C/-)i(A 1 ) J R 21 (A 12 )(i7_) 2 (A 2 ) = (C/_) 2 (A 2 )i?i 2 (A 12 )(C/_) 1 (A 1 )i? 21 (Ai 2 ). 

The corresponding relation for the dual double monodromy matrix is 

J Ri 2 (A 2 i)(^ 1 )i( A i) i? 2i(-Ai 2 - 2i 1 )(U t ?) 2 {\ 2 ) 

= (^ 2 ) 2 (A 2 )i?i 2 (-Ai 2 - 2r ? )([/^) 1 (A 1 ) J R 21 (A 21 ). 

The transfer matrix Txxz(A) € End!K can be constructed from both double row monodromy 
matrices 

Txxz(A) = tr {(K + ) (A) (U-) (X)} = tr {(^) (A)(t^) (A)}. (2.6) 

Here the trace is taken over the auxiliary space Vq. It was shown by Sklyanin [25] that the 
transfer matrices commute for any value of the spectral parameter, 

[Txxz(A),Txxz(m)] =0. 

In the homogeneous limit (£ m = for m = 1 , . . . , N) , the Hamiltonian of the open XXZ spin 
chain with most general boundary fields (12. If) can be obtained as the following derivative of the 
transfer matrix (|2.6|) : 

H = ci-^-T X xz(A)| + const, (2.7) 

where 

ci = — 8sinh(<5) sinh(<5 — rf) sinh(£) sinh(C — rj). 

This procedure permits to construct a commuting family of the conserved charges. However 
to construct the eigenstates in the framework of the algebraic Bethe ansatz one needs a reference 
state |0) which is an eigenstate for the operators A± and D± and annihilated by the operators C±, 

C±(A)|0)=0, A ± (A)|0) = o±|0), D±(\)\Q) = d±|0). 

Evidently, for the spin chains with diagonal boundary terms a ferromagnetic state with all the 
spins up is a state with such properties. But spin chain where the total spin component S z is not 
conserved (like the periodic XYZ spin chain which is related to the elliptic 8-vertex model) do 
not posses a similar reference state. Indeed, to diagonalize the Hamiltonian one needs a gauge 
transformation which maps the underlying vertex model into a SOS model as it was pointed 
out by Faddeev and Takhtadjan |26j . In our case, even if the open XXZ bulk Hamiltonian still 
conserve this U(l) symmetry, the boundary terms break it. Thus our next goal is to construct 
a gauge transformation (vertex-face) which leads to a SOS model with a reflecting end where 
the boundary terms are diagonal. 

3 Vertex-face transformation 

The idea to use a gauge transformation to diagonalize the XXZ Hamiltonian with non-diagonal 
boundary terms was first applied in [3] , then in a more algebraic form the vertex- face transfor- 
mation for this case was studied in [30]. Here we propose a simpler form for this transformation 
and a more explicit construction of the eigenstates. 
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3.1 Dynamical Yang-Baxter algebra 

We introduce the following local gauge transformation 

S(X;e)^S(X;6,u) = e x / 2 ^ 1 1 

This matrix satisfies two important properties. 
• It diagonalizes the boundary matrix ET_(A), 

£_(A) = S-^A;* - (,t)K_(X)S(-X;5- C,r), 
where the diagonal matrix /C_(A) is 

/C_(A) = /C_(A; 5, C) = S-^A; 5 - C, r)i^_(A)5(-A; 5 - C, r) 

\ 



/ sinh(<5 - A) 

jj 

V 



sinh(<5 + A) 






sinh(£ — A) 



(3.1) 



sinh(C + A) / 
It is a vertex-face transformation, namely: 

i?i2(Ai 2 )5i(A i; e)s 2 (x 2 ; e - w \) = s 2 (\ 2 ; e)s 1 (x r , e - v* z 2 )n l2 (x 12 - e), 

where we introduced the dynamical i?-matrix 

/sinh(A + 7]) 




(3.2) 



K{X; 




sinh A sinh(0 — rf) sinh 7] sinh(# — A) 



V 



o 





sinh 9 sinh 9 

sinh 7] sinh(# + A) sinh A sinh(# + rj) 



sinh# 




sinh( 






sinh(A + rf)J 



. (3.3) 



This transformation can be understood as mapping the vertex configuration into a face con- 
figuration in the dual lattice [2] and the i?-matrix fj3.3[) contains the statistical weights of the 
trigonometric SOS model. It is a very well studied object generating the Felder's dynamical 
Yang-Baxter algebra [8]. It solves the dynamical Yang-Baxter equation (DYBE), 

K 12 {X 12 ; 9 - r)a$)K 13 (X 13 ; 6)K 23 (X 23 ; 9 - rjaf) = K 23 (X 23 ; 9)TZ 13 (X 13 ; 9 - ria z 2 )Tl 12 {X 12 ; 9), 

and satisfies the following properties which will be essential to construct the double row mono- 
dromy matrices: 

• Ice rule 

[a z 1 +al,K 12 (X;6)]=0. 







This symmetry is responsible of the six vertex texture of the statistical weight: TV^S 
unless a + f3 = /i + v. It is easy to see that this relation induces a similar relation for the 
transposed i?-matrix: 

[ai-a z 2 ,K t MX;9)]=0, 



which translate into: (T^) 1 ^ = unless a — /3 = fj, — v. 
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• Unitarity: 

K 12 (X; 9)TZ 21 (-X; 9) = - sinh(A - rj) sinh(A + 77) Id. 

• Crossing symmetries. The crossing relations for the dynamical i?-matrix are not as simple 
as for the vertex type i?-matrix, here we write it in the following compact formQ: 

-a\lZ%{-\ -^9 + naOa f^^fH = 7MA; 9), (3.4) 

-a\n^(-X -r,;9- wlH ^^f^ = U 12 (X; 9). (3.5) 

We should stress that these two relations are not equivalent at all since the dynamical 
R- matrix are not symmetric, but there is a weaker condition 

• Parity: 

TZ 2l (X;9) = a x ^a^ v n l2 (\ ] 9)al' y al' y = K 12 (X;-9). (3.6) 

Now we can easily introduce the dynamical monodromy matrix as an ordered product of the 
dynamical .R-matrices 



%(\;9)=n Q1 ^A-£i;0-77X>/j • • -7^jv(A - faO). 



(3.7) 



7o(X;9) satisfies the zero weight condition: 
[T (X;9),a z + S z ] =0. 

It is easy to show that this monodromy matrix satisfies the dynamical Yang-Baxter relation 
K 12 (\ l2 ; 9 - 77S*)Ti(Ai; 9)T 2 (\ 2 ; 9 - w f) = T 2 (A 2 ; 9)T 1 (X l ;9 - W z 2 )TZ l2 (X 12 ; 9), 

N 

where the z component of the total spin S z = ^ af. 

i=l 

This dynamical monodromy matrix is related to the usual one (|2.4p by the vertex-face trans- 
formation S, 

S_({£}; 9)S (X; 9 - V S Z )%(X; 9) = T (X)S (X; 9)S-({£}; 9 - w z ), 

where 



S_({e}; 9) = S N (£ N ; 9) ■ ■ ■ S 1 (^i; 9 - r, ^ a^j 



We need to introduce a second vertex-face transformation to proceed with the construction 
of the double row dynamical monodromy matrix (it will be necessary to construct a dynamical 
analog of the [/+ double row monodromy matrix in the next section). As for the vertex case, 
dual algebras in dynamical context are closely related to the antipode of the corresponding 
dynamical quantum group, which, due to the complicated form of the crossing symmetries, is 
not represented by the simple matrix transposition or inversion. 



1 We assume in this paper the following normal ordering: the ctq in the argument of the 7£oi matrix (which 
does not commute with it) is always on the right of all other operators involved in the definition of 1Z. 
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The matrix 1Z(X; 9) can be also obtained by the second vertex-face transformation defined 
by the same matrix S 

J Ri 2 (Ai 2 )5 2 (A 2 ; 0)5i(Ai; 9 + w z 2 ) = S 1 (X 1 ;8)S 2 (X 2 ; 8 + W l)K 12 {X l2 - 9), 

which leads to the following dynamical Yang-Baxter equation 

n 12 (X 12 ; 8)K 13 (X 13 ; 8 + ^)^ 23 (A 23 ; 8) 

= K 23 (X 23 ; 8 + i7<7f)tt 13 (Ai 3 ; 8)K 12 (X 12 ; 8 + V a z 3 ). (3.8) 

Now we can introduce the "crossed" L-operator, which will be used to construct the dynamical 
analog of the double row monodromy matrix U+ 

and its inverse (in a convenient for us normalization) 

rtu\.ft\- V tu\. ft ^ smh(g + 7?a|) 
L 21 {X,8) - K 21 (A,9 -rjcr^ . 

Due to the dynamical Yang-Baxter equation (j3.8|) and the crossing symmetry (|3.4p . this new 
L-operator also satisfies a dynamical Yang-Baxter type equation: 

Tl 12 (X 2l ;8 + W Z 3 )C[\(X 13 ; 8)£ 2 %(X 23 ; 8 - W z ) = C%{\ 23 ; 9)C%{X l3 - 8 - ryaf )fti 2 (A 21 ; 8). 

It is easy to derive the properties of this L-operator, such as unitarity 

£ 2 \(-A - m 8)C%{\ - V ;9) = - sinh(A - ri) sinh(A + rj) Id, 

ice rule 

[af - a z 2 ,C%(X;9)} = [of - a z 2 ,^(X;9)} = 0, 
and parity 

af' y a 2 ' y C^ 2 (X; 9)a^ y a 2 ,y = C 2 \(X;9) = £%(X; -8), 

from the corresponding symmetries of the dynamical i?-matrix. 

Starting from this L-operator we can construct a dual monodromy matrix 



V*° (A; 9) = £*V A - &v; 9 + rj £ af - ■ ■ £*(A - £ i; 9). (3.9) 




The Yang-Baxter relation for this matrix can be written in the following form 

n 12 (X 2l ;8 + r/S^Vj 1 (Ai; 0)V£(A 2 ; 8 - W \) = V* 2 (A 2 ; 9)V\ 1 (X 1 ; 9 - r 1 a z )n 12 {X 21 ; 9). 

Once again, it is possible to construct a gauge transformation which relates this matrix to the 
transposed vertex monodromy matrix 

£+({£}; 9)S(X + 77; 8 + ^) V*° (A; 8) = T*%X)S(X + 77; 8) S + ({C}; 8 - w z ), 

where 

S + m; 0) = Si(&; 9) ■ ■ ■ S N Un; O + V^fj, 

and 

S (X;9) = a y S (X;9)a y . 
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3.2 Dynamical reflection algebra 

Now we can proceed with a construction of the double row monodromy matrix. The key ingre- 
dient is once again the reflection equation. It easy to check that for the diagonal /C matrix fj3. If) 
satisfy the following dynamical reflection relation for 9 = 5 — C, 

fti 2 (Ai 2 ; 0)(£-)i(Ai)^2i(Ai 2 ; 0)(/C-) 2 (A 2 ) 

= (K-) 2 (A 2 )ft 12 (A 12 ; 0)(£_)i(Ai)ft 2 i(Ai 2 ; 0). (3.10) 

This is essentially the reflection equation introduced in [4], with the dynamical i?-matrix. This 
equation can be checked directly but also it can be derived from the vertex reflection equation 
using the gauge transformation S: 

Proposition 3.1. The vertex type reflection equation (|2.5p for K- is equivalent to the SOS type 
reflection equation for JC- if we the SOS parameter is chosen as9 = 5 — C,,uj = t. 

Using (|3.2p the left hand side of the reflection equation (|2.5p can be rewritten as 

i?i 2 (Ai 2 )if_(Ai)i? 2 i(Ai 2 )ir_(A 2 ) = s 2 (x 2] 9)s l (x l -9-r i a z 2 )n l2 {x l2 ;9)ic l (x l ) 

x K 2l {X l2 - 6)K 2 (\ 2 )S^{-\i;e - W Z 2 )S^(-X 2 ; 9). 

To obtain this formula we used the fact that /C_ is diagonal. Similar calculation for the right 
hand side give us: 

K-(A 2 )E 12 (Ai 2 )if_(A 1 ) J R 21 (A 12 ) = 5 2 (A 2 ;^)5 1 (Ai;0-r ? (rf) JC 2 (X 2 )n 12 (X l2 ; 9) 
x /C x (A i; )K 21 (X 12 ; 9)S^ 1 (-\ 1 ;9 - W Z 2 )S^(-X 2 ; 9). 

Which finally leads to the equation (|3.10p . In a similar way one can prove that the dynamical 
double row monodromy matrix defined as 

(U.) (X;9) = r o (A;0)(/C_) o (A)fo(A;0) = (^^9) pI(A;S)) „ ' (3 - U) 

where 

= 7(A)og7?»(-A - m + ^)a f nh ^p = tWV^-A; 9), 

satisfies the dynamical reflection algebra relation, 

^i 2 (Ai 2 ; 9 - n S*)(W_)i(Ai; 9)K 21 {X 12 9 - r,S z )(U-) 2 (X 2 ; 9) 

= (W_) 2 (A 2 ; 9)TZ 12 (X l2 ; 9 - r, S 2 )(W_)i(Ai; 9)TZ 2l (X 12 ; 9 - r, S 2 ). (3.12) 

This double row monodromy matrix has a clear interpretation in statistical mechanics, it is 
an essential tool to study the trigonometric SOS model with a reflecting end. We will return to 
this correspondence in the last section of this paper. 

It can be easily shown that the vertex-face transformation S gives a relation between this 
dynamical double row monodromy matrix and the usual one 

£_({£}; 9)S (X; 9 - V S Z )(U^) (X; 9) = (£/_)o(A)S_({£}; 6)S Q (-\; 9 - rjS z ). (3.13) 
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3.3 Dual reflection algebra 

In the vertex case we constructed two double row monodromy matrices: one around the K_ 
matrix and second one around the K + matrix. The same procedure is possible in the dynamical 
case. The only difference is that because of a more complicated crossing relations (|3.4p and (j3.5[) 
instead of the transposed monodromy matrix we will use here the dual monodromy matrix 

v*°(a) (prgp . 

We start once again with a diagonal solution /C+(A) of the dual reflection equation with the 
dynamical i?-matrix, 

n l2 {\ 2l -'6){lC%) l {\ l )n 21 {-X l2 - 2?7;£)(/C' + 2 ) 2 (A 2 ) 

= (/C i + 2 ) 2 (A 2 )^ 12 (-Ai 2 - 2 V ;9)(}C t l) 1 (X 1 )n 2l (\ 2l - } 9). 

It can be expressed in terms of the fC- matrix, 

/C + (A) = /C + (W) = /C_(-A - V ;S,C). 

The crucial point here is that both the vertex boundary matrix K + and the dynamical bounda- 
ry /C+ are related each other by mean of the second Vertex-face transformation: 

£*_(A) = S~\\ + m S- C, r)if*.(A)5(-A - m S- C,r). 

Note that the dynamical parameter 9 should not a priori coincide with the parameter 9 from 
the previous section as it is chosen here 6 = 5 — £ to diagonalize the boundary matrix K + . 
As in the previous section we can construct the following double row monodromy matrix 

( M+)S '(A;S) = VS(A;5)(^)„(A)*(A;S) = ^ . 

where 

V*°(A; 0) = £&(A + &;0) ■ ■ ■ £% (x + Cme + ^J^ < ] = i(\)(V$>y\-\ - 2 m 9). 



i=l 



The normalization factor in the last equation can be easily computed using the unitarity of the 
operator £, 



v 



7(A) = 7(A) J] sinh(A + & ) sinh(A + & + 2rj). 



i=l 

This double row monodromy matrix satisfies the following dynamical reflection equation 

Tl 12 (X 21 ; 9 + r?S z )(^ )j (Ai ; 0)^21 ( - A i2 - 277; 5 + 17 S 2 )) (W*? ) 2 (A 2 ; 0) 

= K 2 ) 2 ( A 2 ; 5)7^12 (- Au - 2 ??; + r, S Z )(U^ )i (Ai; 0)ft 2 i (A 2i ; £ + 77 S z ) . (3. 14) 

Its relation with the corresponding vertex double row monodromy matrix can be easily written 
in the following form 

S + aCy,0)S (X + 7 1 -,9 + i ] S z ){U t + ) (X-,9) = (U t + ) Q (X)S + aO;0)S (-X-r ] ;9 + r ] S z ). 

In the next section we show how to use these matrices to construct eigenstates of the open 
XXZ spin chain with non-diagonal boundary terms. 
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4 Algebraic Bethe ansatz 

Our goal is to construct eigenstates of the spin chain Hamiltonian (|2.1|) which is equivalent 
to the construction of eigenstates of the transfer matrix Txxz(A) ()2.6|) . We have seen that 
usual algebraic Bethe ansatz procedure cannot be applied to this case as there is no simple 
reference state. Thus our strategy is to consider dynamical monodromy matrix instead of vertex 
one using the vertex-face correspondence. The advantage of this approach is the fact that this 
procedure allows us (under certain constraints) to diagonalize the boundary matrices, and hence, 
it becomes possible to use the algebraic Bethe ansatz techniques. 

Using (|3.13j) the transfer matrix (|2.6p can rewritten in terms of the dynamical algebra gene- 
rators 



This trace can be rewritten in terms of the operator entries of the double row monodromy matrix 



Tr {K + (A)5 (A; 9 - V S Z )(T-) (X; 0)Sq~ (— A; 9 - r,S z )} 
= {So\-X; 6 - r)S z )K + (X)S (X; 9 - r)S z )}tA-(X; 9) 
+ {S^i-X; 9 - V S Z )K + (X)S (X; 9 - V S Z )}ZV„(X; 9) 
+ {SoH-X; 9- v S z - 2r,)K+(X)S (X; 9 - r)S z )}+B-(X; 9) 
+ {So\-X; 9-r]S z + 2 V )K + (X)S (X; 9 - r]S z )}±C-(X; 9). 



The next crucial step is to restrict our analysis to a subspace with a fixed z component of the 
total spin (in the face picture, it has nothing to do with the XXZ spin). This restriction follows 
from the results of [3[ [21] , it was shown that the number of Bethe roots for all the eigenstates is 
the same. It means that in the framework of the algebraic Bethe ansatz the number of creation 
operators and hence the total spin should be fixed. Unfortunately it also means that in this 
approach the algebraic Bethe ansatz cannot lead to the complete description of the eigenstates. 

Thus we consider the action of the trace in (|4.ip on the states with a fixed z component 
of the total spin 



To apply the algebraic Bethe ansatz we need to diagonalize both boundary matrices simul- 
taneously, it means that the non diagonal terms should be zero, 



T xxz (X)S-m,9) = S-({Z},9) 

x Tt {K + (X)S (X; 9 - V S')(U-) (\; 9)S« 1 (-X; 9 - V S Z )}. 



(4.1) 



s 2 IV> 



{So '(-A; S-C-m, r)K + (X)So(X; 5 - ( - r,(s - 2), r)} + = 0, 
{Sq-^-A; 5 - ( - vs,t)K + (X)S (X; 8-(- V (s + 2),r)}+ = 0. 



These lead to two constraints on the boundary parameters: 



cosh(5 — C) = cosh((5 — C — r]s + t — r — rf) 
cosh(<5 — C) = cosh(<5 — £ — — r + r + rf) 



(4.2) 
(4.3) 



which can be solved by imposing: 



r = t + 7] + ira 



5 — ( = 5 — £ — rjs + 2mm. 



(4.4) 



With these constraints we obtain the following result: 



Txxz(A)S„({£}; 0)W) = S_({£}; 0)Tr o (/C + (A; 6, C)71(A; 0)M, 
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where 

smn(o — (J 

Therefore the eigenstates of the vertex transfer matrix Txxz can be obtained from the eigen- 
states of the SOS transfer matrix 

T SOSl (A; 9) = Tr (K+ (A; 6, QU- (X;9)). (4.5) 

Using the reflection algebra relation (]3. 12|) which contains the commutation relations for 
the operators A-(X;9), B-(X;9), C-(X;9), and T>_(\;9), we can actually construct two sets of 
eigenstates of TsoSi (A; 9). 

The advantage of the SOS picture is the existence of the reference state |0) with all the spins 
up, which is an eigenstate of the operators A-(X;9) and T>_(X;9) and it is annihilated by the 
operator C-(X;9), C-(X;9)\0) = 0. Then the eigenstates are constructed using the algebraic 
Bethe ansatz technique: 

Proposition 4.1. Let 9 = 5 — Q, and 5, Q, 5 and £ satisfy the boundary constraints (|4.4|) with 
total spin s being even if N is even and odd if N is odd, \s\ < N. Then the state 

N — s 

\t/> 1 -({X}m))=B-(\ 1 ;0) — B-(\ m ;9)\0) 1 M = — ^— , (4.6) 

is an eigenstate of TgoSi 9) with eigenvalue 

sinh(£ — fi) sinh(<5 + fi) sinh(<5 — /i) sinh(2^ + 2rj) 



Ax (n; {A}; 6, (, 5,() 



sinh(C — /j, — rj) sinh(<5 — fx — rj) sinh(c) + fi) sinh(2 / u + rj) 

n 1 sinh(/i + Xi) sinh(u — A; — r?) -r-r , . . . . , . 



+ 



. J- sinh(ii + Xi + rj) sinh(/x - Aj . 

sinh(£ + /i + rf) sinh(<5 + jjl + rj) sinh(£ — fx — rj) sinh 2fi 
sinh(C — [i — rj) sinh(<5 + fx) sinh(£ + fx) sinh(2^ + rj) 

Jt- sinh(/i + Xi + 2rj) sinh(/i — Xj + rj) -^ lM ■ u c n 

x 11 — r r? — r~\ — i — r~~T7 — 11 smh (A* + 6 smh )• 

"7 smh( / u + Aj + r/) smh(/x - Ai) A ^ 
i=i i=i 

provided the Xi, i = 1, . . . , M satisfy the Bethe equation: 

y 1 (X i ,{Xy,6,(,6,0 = y 1 (-X i -ri,{X};6,(,5,0, i = l,...,M, (4.7) 

yi(Aj, {A}; J, £5, C) = sinh(C + Aj)sinh(5 - Aj)sinh(C - Aj) sinh(<5 + Aj) 

71/ /V 

x sinh(Aj + X k ) sinh(Aj - A fc - rj) x sinh(A.j + Cj + fj) sinh(Aj - ^ + rj). 

k^i j=l 

The proof of this proposition is a standard algebraic Bethe ansatz computation, some details 
on the commutation relations leading to this result are given in the Appendix [Aj 

It is possible to construct a second set of Bethe states starting from the second reference 
state |0) with all the spins down. It is also an eigenstate of the operators A-(X; 9) and T>_(X; 8), 
while S_(A;0)|O) = 0. 
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Proposition 4.2. 

IV- ({\}m)) = C_(Ai; 6) ■ ■ ■ C_(A A/ ; 9)\% M = (4.8) 
is an eigenstate o/Tgosi(/^;0) with eigenvalue: 

A 2 (/i;{A};J,CJ,C) = Ai(A t ;{A};C,<5,CJ), (4.9) 
provided the Ai, i = 1, . . . , M satisfy the Bethe equations: 

y 2 (\ i ,{\};5,(,8,0 = y 2 (-\ i - V ,{\};5,(,5,(), i = l,...,M, (4.10) 

with 

Ite(Ai,{A};<5,Cj,C) = yi(Ai,{A};C,*,C,*)- (4.11) 

These two set of Bethe states are actually related each other. More details on this relations 
are given the Appendix [Bj Now it is easy to construct the eigenstates of the XXZ transfer 
matrix; 

Proposition 4.3. Let 9 = 5 — Q, and 5, 5 and £ satisfy the boundary constraints (|4.4j) with 
total spin s being even if N is even and odd if N is odd, \s\ < N. Then the states 

S_(tf};Wi({A}M B )>, 

with Mb = ^y^- are eigenstates of Txxz(A) provided \i,...,\m b satisfy the Bethe equa- 
tions KTb . 

A second set of Bethe states can be constructed as 
S_({£};0)h^({A}M c )>, 
with M c = if Ai, . . . , A A / C satisfy (|4TUj) . 

The corresponding eigenvalues of the Hamiltonian is then given by (}2T 

M 



E = y^e{\j) + ciA^coth?], e(Aj) 



ci sinh ?] 



sinh(A + r/) sinh A 

Thus we obtained a simple algebraic expression for the eigenstates of the spin chain with non- 
diagonal boundaries in terms of the SOS operators and vertex-face transformation S. This 
method however has two clear limitations: first, it works in its present form only if the con- 
straints (|4.4p are satisfied, and second, it permits to construct only states with fixed number of 
operators £>_(A;#) (or C_(A;0)) and this set of states is evidently incomplete (it is clear from 
the SOS picture). Thus the problem of the complete description of the spectrum remains open. 

The construction of the same eigenstates can be equivalently performed starting from the 
U+(\;0) reflection algebra. Following the same lines as before we obtain: 

T xxz (x)s + m-,e) = s+m-,9) 

xTro{5o(A + 7 ? ;5 + r ? S z )(W| o )o(A;0)^ 1 (-A-7 7 ;e + 7 ? S^i(A)}. 



Again, it is possible to write the trace in terms of the operator entries of the double row monod- 
romy matrix the W + (A, &). It is easy to check that if we consider the action of this trace on the 
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states with z component of total spin s and if the constraints (|4.4p are satisfied the non-diagonal 
terms in this expressions vanish 

Txxz(A)S + ({£}; = £+({£}; #)Tr (^ (A; 0)£*_°(A; 5, ()) 

where 

? A jf ^ sinh((5 - C - W z ) y / \ x ^ 

£_(A; 5, C) = ; 7—77 — £+(-A - 77; 5, C). 

smn(o — Q 

It means that we replace once again the vertex transfer matrix with non-diagonal boundary 
matrices by a SOS transfer matrix with diagonal boundaries. 

T S os 2 (A, 9) = Tr («f (A; 9)1?° (A; 5, C)) • 

Proposition 4.4. Lei = 5 — and 5, 5 and £ satisfy the boundary constraints (|4.4p wt/t 
ioiaZ spin s freing even i/iV is e?;en and odd if N is odd, \s\ < N. 

N — s 

IV>!({A}m)> = 6+(Ai; 0) • • • S+(A M ; 6>)|0>, M = — ^— , 

is an eigenstate 0/ Tsos 2 (/•*; ^i/i eigenvalue A2(//; 0", £, <5, Q provided the \, i = 1,...,M 
satisfy the Bethe equations 

y 2 (A i ,{A};5,C,5,C) = y2(-A i -n,{A};5,C, ( 5,C), i = l,...,M. (4.12) 
Proposition 4.5. 

|^({A})) = c + (Ai;0) • • -c + (A M ; 0)|o), m = ^±i, 

is an eigenstate o/Tsos 2 (w#) eigenvalue Ai(/x; o~, £, 5, Q provided the Aj ; i = 1, ...,M 

satisfy the Bethe equations 

yi(Ai,{A};I,C,*,C)=yi(-Ai-»7,{A};5,C,«y,C), i = l,...,M. (4.13) 

Finally we can once again construct the corresponding eigenstates for the XXZ transfer matrix 
using the vertex- face transformation S+. 

Proposition 4.6. Let = 5 — £, and 5, £ ; <5 and £ satisfy the boundary constraints (|4.4p uiii/i 
toto/ spin s freing even i/iV is e?;en and odd i/Af is odd, \s\ < N. Then the states 

s + m-,9)\^ + ({x} MB )), 

with Mb = are eigenstates of Txxz(A) provided \i,...,\m b satisfy the Bethe equa- 

tions (gig) . 

^4 second set of Bethe states can be constructed as 

s+({£};W+({a}m )>, 

with M c = ^y 1 , if Ai, . . . , Am c safc/y (|41~3jl . 
It is easy to check that 

Al(n;5,CJ,C) = A 2 (/i;5,C,«5,C), 
y 1 (A i ,{A};5,C,5,C) = n 2 (A 4 ,{A};5,C, ( 5,C). 

It means that the eigenstates constructed by the first method (Proposition I4.3P and the states 
constructed using the dual algebra (Proposition I4.6P are exactly the same states. However it was 
already shown for diagonal boundary conditions [16] that it is important to be able to construct 
the same states by two different methods to proceed with the computation of the scalar products 
and correlation functions. 
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5 SOS model with reflecting end 

In the previous section we used the trigonometric dynamical reflection algebra to construct the 
eigenstates of the spin chains with non-diagonal boundary fields. However this algebra has 
another clear interpretation, it describes a trigonometric solid-on-solid (SOS) model with one 
reflecting end (in the same way as usual reflection algebra with diagonal matrix K describes 
a six- vertex model with reflecting end). 

The SOS model is a two dimensional statistical mechanics lattice model which can be defined 
in terms of a height function. Every square of the lattice is characterized by a height 9 and its 
values for two adjacent squares differ by i). There are 6 possible face configurations 



9-7] 


9-2n 


9 + T] 


9 + 2r] 


9-7] 


9 


9 


9-7] 


9 


9 + T] 


9 


9 + 7] 


9 + 7] 


9 


9 + 7] 


9 


9-7] 


9 


9 


9-n 9 


9 + t] 9 


6-7) 



and the corresponding statistical weights are given by the dynamical iZ-matrix (13.3p . 

We consider this model with a reflecting end, which means that each horizontal line makes 
a U-turn on the left (right) side of the lattice. It produces two following configurations charac- 
terized by the boundary matrix /C_(A) (fC + (\)): 




IC-(X)u /C_(A) 22 
Figure 1. Boundary configuration with external height 6. 

It was seen in the vertex case that it is crucial to study the partition function of this model 
with domain wall boundary conditions. It is, in particular, necessary to study the correlation 
functions of the spin chains with non-diagonal boundary terms. 

These domain wall boundary conditions can be easily written for both cases: with left or right 
reflection end, here we give the example of the left reflection end (given by the matrix /C_), see 
Fig. [2 Here the heights decrease from left to rights on the upper boundary, the heights grow from 
left to right on the lower boundary (we will call this situation case I). As left external height 
is fixed these two conditions determine completely the configuration on the right boundary 
(heights decreasing in the upward direction). Equivalently one can consider heights increasing 
on the upper boundary and decreasing on the lower one (case II). The same two situations (case I 
and case II) are possible for the reflection on the right boundary (with boundary weights given 
by the matrix 

Thus we have four possible configurations and for all of them the partition functions can be 
written in terms of the double row monodromy matrix. Namely for the left reflecting end the 
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8-Tl 




9-(N-l)ti 


0-Nii 








8-(N-1)ti 








e-(N-2)r| 








6-(N-3)r| 






















e+(N-3)ii 








9+(N-2)r| 








6+(n-i)ti 


e+ri 




9+(N-l)ti 


8+N11 



Figure 2. Domain wall boundary conditions. 



case I partition function can be written as 

N 

4; 2 iv({A}, {£}, 5, C) = <o| n S -(A*; 9)\0), 

i=l 

while for the case II we obtain 

N 

i=i 

where we used our standard notations 6 = 5 — £. 

If the reflecting end is on the right, we obtain the following representations: 

N 
N 

Z C N + 2N ({\}, = (0|nc+(Ai;5)|0>. 

i=l 

All these partition functions can be computed following the same lines. The detailed derivation 
of the function Z N ~ 2N is given in [TTj , it is shown that this function is uniquely defined by the 
following properties 

• For each parameter Aj the normalized partition function 

Z n ,2n({\}, {£}, 6, = exp ($N + 2) 

x sinh(5 + A<) sinh(C + X l )Z^ 2N ({\}, {£}, S, (), 
is a polynomial of degree at most 2N + 2 in e 2Ai . 
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For N = 1 the partition function is 



ryts / \ j- e >\ sinh7?sinh(5 - ( - r\) /sinh(5 - A) 

zf 2 (a, e, 5, o = — ^-727^ — ^ — - ( „^) s , x smh ( A - o smh (^ - c + a + e) 



+ 



sinh 2 ((5-C) V sinh (^ + A ) 
sinh(C — A) 



sinh(C + A) 



sinh(A + £) sinh(<5 - C ~ A + £) 



• 2 jv,2iv({A}> C) is symmetric in A*. 

• ^a^at^A}, C) is symmetric in 

• Crossing symmetry 



Z N,2N(~ X i ~ V, {A}, {£}> 
sinh(Aj + 5) « 



sinh(2(Ai + t/)) sinh(A; + C) 



sinh(Aj — 5 + 77) 



• Recursive relations 



sinh(2Aj) sinh(Aj — ( + r/) 

^(Ai,{A},{e},«,0- 

sinhrysinh((^ — Ai) 
A 1= 5i sinh(C + Ai) 



TT • ^< , smh(6 - ( + (N - 2i) V ) 
n Smh(A ' +6) sinh(^-C + (^-2, + l),) 



iV 



x sinh(Ai - & + 77) sinh(Ai + & + 77) sinh(A.; - £1 + 77) 

i=2 

X Z (N-1),2(N-1) ({A}2,...,JV, {£}2,...,AT> ^ C)> 
4;27V({A},{£WC) 



sinh 77 sinh((5 — An) 



A? 



JsinhCAi-a) 



Ajv=-£i sinh(5 + A^) 

sinh(<5 - C + (JV - 2i)rj) 



sinh(<5 - C + (N - 2i + 1)??) 
x JJ sinh(AAr + & + 77) sinh(AAr - £i + rj) sinh(Ai_i + £1 + 77) 

i=2 

x ^(JV-i),2(AT-l)({A}l,...,Ar-l, {€h,...,N,S,0- 



i=l 

N 



It means that if we find a function satisfying all these properties it is the partition function. 
The following proposition can be easily proved by the direct verification of these properties: 

Proposition 5.1. 



v 



z n~2n({ x }> = 7(A) det M i:i J 



sinh(<5 - C + rj(N - 2i) 
L± \smh{5 - C + r]{N - i)) 



N 



n sinh(Aj + £j) sinh(Aj - £,) sinh(Aj + £, + 77) sinh(Aj - £j + 77) 

i,j=l 

sinh(£j + £j) sinh(£j — £j) sinh(Aj — Aj) sinh(Aj + Aj + 77) 

l<i<j<N 
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where 

= sinh(<5 + g,,) sinh(C-gj) 
i,j sinh(<5 + Xi) sinh(C + A*) 



sinh(2Aj) sinh rj 



sinh(Aj — £j + rj) sinh(Aj + £j + rj) sinh(Aj — £j) sinh(Aj + 
The partition function Z^~ 2N can be computed in a similar way: 

Z C N - 2N ({X},m,5,C) = Z B N - 2N ({X},{^,C,S). (5.1) 

It is also possible to derive it from the parity symmetry (|3,6p of the dynamical i?-matrix. More 
details on this are given in the Appendix 151 

The other two partition functions can be easily obtained using the same method. It is 
also possible to derive them from the two previous results using an isomorphism between the 
dynamical reflection algebras generated by the elements of IA~ and U+. We give more detail on 
this isomorphism in the Appendix O Using this isomorphism we can easily prove the following 
propositions: 

Proposition 5.2. 

4 + «A}, {£} J,C) = z%-{{ -x - V }, 

Proposition 5.3. 

Z B N + ({A}, {£}, 5, = Z C N ({ -A - rj), {£},!, (). 
Notice also that: 

4 + ({A},{e}J,o = 4 + «A},{e},cJ). 

We would like to mention that similar results can be obtained from the so-called .F-basis repre- 
sentation for the B± and C± operators. This idea was first developed in [18] for the six-vertex 
model. The F-basis for the SOS model was constructed in [1] and its application to the case 
with a reflecting end was proposed in [31J . The computation of the partition function using this 
approach was recently presented by W.L. Yang [29] , 

These results are interesting from two different points of view. First, it is the first case where 
the partition function of a SOS model with domain wall boundary conditions is computed as 
a single determinant (and not as a sum of determinants as in [21]). Second, these expressions, 
as it was shown in the previous sections, can be directly used to study the XXZ model with 
non-diagonal boundary terms. In particular, it is the first step towards the computation of the 
correlation functions and form- factors. 

It is important to underline in conclusion that we considered here a very special case of the 
spin chain with non-diagonal boundaries: the boundary magnetic fields are not independent 
but satisfy two boundary constraints (|4.4p . This case is interesting as it leads to a direct 
correspondence with a SOS model with diagonal boundary matrices. However from the spin 
chain point of view it is a degenerated situation which is in some sense more complicated 
to study than the general case with only one constraint (either condition (14. 2\i or (14. 3p ). In 
particular, the completeness of the Bethe ansatz conjectured for the latter case [22] is evidently 
missing in the situation considered in this paper. Thus, there are two important open problems: 
construction of the missing states in the case with two constraints and application of the present 
approach to the most general case solvable by Bethe ansatz with only one constraint. 
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A Dynamical commutation relations 

In this appendix we give the derivation of the dynamical algebraic Bethe ansatz. We prove here 
a more general version of Proposition 14. 11 Propositions 14.2 1 14.41 and 14.51 can be proved in a very 
similar way. 

We consider here the following dynamical transfer matrix 

T SOSl (A; 0) = Tr (£+ (A; 9 + ( - nS z , QU- (A; 9)) . 

After our usual restriction on the subspace with the fixed z component of the total spin and 
using the boundary constraints (I4.4p we obtain exactly the same transfer matrix (14.50 as in 
Proposition 14.11 

Now we can construct the eigenstates of this transfer matrix, following the usual boundary 
algebraic Bethe ansatz technique [25J. 

First it is convenient to introduce a modified operator Z?_(A;0), 

~ sinh(6l - rjS z + rj) ( _ sinh(6l - riS z + 2A + rj) sinh 77 

1 ' j sinh^-^) \ U -W> sinh(2A + 7 ? )smh(0-r ? S z + r ? ) 1 ' 

It is easy to compute the action of the operators A-(X; 9), £>_(A; 0) on the reference state |0) 

A- (A; f ) |0) = Sm ^~^ J] sinh ( A + sinh ( A +€i + »?)|0>, 
sinn(o + A) x - 1 - 



i=l 

N 



^ /> sinh 2A sinh(C — A — ??) sinh(£ + A + ??) tt . , ^ . 

P - A; = • uo^x \ ■ u/lu ■ ux^ vT II smh A - & smh A + f ° • 
smh(2A + 77) sinli(( + A) smh(d + A) 

The dynamical reflection relation (I3.12p gives the following commutations rules for the opera- 
tors A-, T>- and B- 

sinn(ft — 77a 2 — 77) smn(Ai 2 + 77) 

smh(Ai 2 ) smh(Ai 2 + 77) 
sinh, sinh 2A 2 sinh(A 12 - 9 + + v) B _ (Ai; ^ (Aa;(?)j 



P_(A i; 0)£_(A 2 ;, 



sinh(6> — 77S Z — 77) sinh(Ai 2 ) sinh(2A 2 + 7/) 
sinh(Ai 2 + 9 - rjS z ) 
sinh(# — 77 S z — rj) 

sinh 77 sinh 2A 2 sinh(2Ai + 277) 



B_(Ai;0).A_(A 2 ;i 



/3_(A i; 0)P_(A 2 ;0). 



sinh(Ai 2 + 77) sinh(2Ai + 77) sinh(2A 2 + 77) 
+ sinh(A 12 + r?) sinh(A 12 + 277) ~ 
sinh(Ai 2 ) sinh(Ai 2 + 77) 
sinh 77 sinh(2(Ai + 77)) sinh(Ai 2 + 9 — t]S z — if) 
sinh(Ai 2 ) sinh(2Ai + 77) sinh(# — 77s 2 — 77) 

The transfer matrix can be expressed as 

T SOSl (A;S> = i?§±A±^_ ( A;*) 
sinn(( — A — 777 

sinh(C - A) sinh(C + 9 - rjS z + A) sinh(2A + 277) 

H = = ^4_(A; 9). 

sinh(£ — A — 77) sinh(£ + 9 — rjS z — A — 77) sinh(2A + 77) 
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Now easily one can check using usual algebraic Bethe ansatz arguments that a state constructed 
by the action of operators £>_ 

\il>l({\}M))=B-(\i;e)---B-(\ M ;e)\0), 

for < M < N is an eigenstate of the transfer matrix TsoSiO- 4 )^) provided the spectral 
parameters satisfy the Bethe equations (I4.7|) with 5 = 5 — C + C — v(N ~ 2M). Note that 
this result is slightly more general than Proposition I4.lt as we can construct eigenstates with 
arbitrary value of the z component of the total spin iV — 1M. However it is important to 
remember that only the eigenstates with a fixed value of M are related to the eigenstates of the 
spin chains with non-diagonal boundary terms. 

B Relations between Bethe states 

A simple relation exist between the two set of Bethe states (I4.6D and (14. 8D by mean of a gene- 
ralized parity symmetry of the dynamical double row monodromy matrix (|3.1ip . With the help 
of the parity symmetries ()3.6p for the dynamical i?-matrix, we easily find the corresponding 
symmetry for the dynamical monodromy matrix ()3.T() : 

o%T(\\6)o% = T X T(\;-6)T X , 

N 

where T x = Y\crf,we use once again our notation 9 = 5 — £. Namely, parity symmetry is 

i=l 

equivalent to exchanging 5 and £. 

A similar relation exist for the solution ()3. 1 [> : 

^/C_(A;J,CK = /C_(A;C,<5), 
leading to the parity symmetry for U~(\; 5 — Q: 

a x U-(\;5- CK = T X U.(X;C- S)T X , (B.l) 
which implies the relations between B_(\, —6) and C_(A,#): 

C-(X,S-0 = T X B-(X,C-5)T X . 
Finally, using (jB.ip . it is obvious that: 

T SO s 1 (A;e) = Tsosi(A;*,CJ,C) = rxT S osi(A;C,<y,C,*)r a! . (B.2) 

This last relation (|ET2jl gives a clear understanding of the relations (|4.9p and ([4. lip between 
the two set of Bethe states, their eigenvalues and Bethe equations. Namely, the two set are 
actually related by interchanging 5 and £. This hidden parity symmetry from the SOS point of 
view (9 — > —9) remains obvious if we consider the Hamiltonian (I2.ip of the spin chain, which is 
symmetric in (<$,£) and (5,()- Using this relation, we obtain a simple algebraic derivation of the 
relation ()5.ip for the partitions functions Zj7 2 n an d %n~2N- 

C Isomorphism between the dynamical algebras 

We define the following isomorphism between algebras (|3. 12[) and (|3.14p : 



p: T-(\; 9) ^T x 7l{- \- m 0)T. 
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It is easy to demonstrate that p(T~(—X — r];6)) satisfies the algebraic relations (|3.14p . This 
isomorphism tell us that, there is a direct relations between the two a priori different SOS 
models with reflecting end described by f|3. 12|) and (|3.14p . We can use this correspondence to 
establish a simple relation between partitions functions. It is easy to notice that C+(A;#) = 
p(B-(X;9)) =T x B-(-X-v,0)T x , hence 

N N 

(o| Hc + (\ i; d)\o) = <o| JJiVM-Ai - mOWxfi) 

i=l i=l 
N 

<0| n*M-Ai - V,0)\0) = Z B N ({-\ - r)}, {£}J,0- 

i=l 

Similar relations can be obtained for two other partitions functions. This isomorphism gives 
also a clear understanding of the relations between the Bethe states for the XXZ spin chain with 
non diagonal boundaries constructed using two algebras. 
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